Theory of spectrum in qubit-Oscillator systems in the ultrastrong coupHng regime 
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Recent measurement on an LC resonator magnetically coupled to a superconducting 
qubit arXiv:1005.1559 shows that the system operates in the ultra-strong coupling regime and 
crosses the limit of validity for the rotating-wave approximation of the Jaynes-Cummings model. 
By using extended bosonic coherent states, we solve the Jaynes-Cummings model exactly without 
the rotating-wave approximation. Our numerically exact results for the spectrum of the flux qubit 
coupled to the LC resonator are fully consistent with the experimental observations. The smallest 
Bloch-Siegert shift obtained is consistent with that observed in this experiment. In addition, the 
Bloch-Siegert shifts in arbitrary level transitions and for arbitrary coupling constants are predicted. 



PACS numbers: 42.50.Pq, 03.65.Ge, 85.25.Cp, 03.67.Lx 

I. INTRODUCTION 

The Jaynes-Cummings (JC) model[l| describes the in- 
teraction of a two-level atom with a single bosonic mode, 
which is fundamental model in quantum optics. Recently, 
the JC model is also closely related to condensed mat- 
ter physics. It can be realized in some solid-state sys- 
tems, such as one Josephson charge qubit coupling to an 
electromagnetic resonator [2l , the superconducting quan- 
tum interference device coupled with a nanomechanical 
resonator Q, and the LC resonator magnetically coupled 
to a superconducting qubit [S]. In conventional quan- 
tum optics, the coupling between the "natural" two- level 
atom and the single bosonic mode is quite weak, the 
rotating-wave approximation (RWA) has been usually 
employed. With the advent of circuit quantum electro- 
dynamics (QED), on-chip superconducting qubits (the 
"artificial " two- level atoms) could be engineered to in- 
teract very strongly with oscillators (cavities) 2 9] , RWA 
can not describe well the strong coupling regimeflO], so 
the studies to the JC model without RWA is highly called 
for. 

However, it is more difficult to solve the JC model 
without RWA than with RWA. In the absence of RWA, 
due to the presence of the counter-rotating terms, the 
photonic number is not conserved, so the photonic Fock 
space has infinite dimensions. The standard diagonaliza- 
tion procedure (see, for example, Ref. [5j]) is the first can- 
didate, which is to apply a truncation procedure consid- 
ering only a truncated number of photons. Typically, the 
convergence is assumed to be achieved if the numerical 
results are determined within very small relative errors. 
Within this method, one has to diagonalize very large, 
sparse Hamiltonian in strong coupling regime. Further- 
more, the calculation might become prohibitive for higher 
excited states where more photons should be involved. 

Fortunately, several non-RWA approaches |10l - [l5| has 
been recently proposed in a few contexts. Especially, 



by using extended bosonic coherent states, three of the 
present authors and a collaborator have solved the Dicke 
model without RWA exactly in the numerical sensejll|. 
The JC model is just special Dicke model with only one 
two-level atom. 

Recently, the spectrum for an LC resonator magneti- 
cally coupled to a superconducting qubit was measured 
experimentally. A 50 MHz Bloch-Siegert shift when the 
qubit is in its symmetry point was observed, which clearly 
shows that the system enter the ultra-strong coupling 
regime. Therefore JC model with RWA is invalid to de- 
scribe this strong coupling system. In this paper, we nu- 
merically solve JC model without RWA exactly. Based 
on the some key data drew from the spectrum'5|, we ob- 
tain a fit of the experimental parameters. All spectrum 
line can then be calculated. The Bloch-Siegert shifts in 
arbitrary level transitions and in a wide range of the cou- 
pling parameters can also be estimated. 

The paper is organized as follows. In Sec. II, the nu- 
merically exact solution to the JC model is proposed in 
detail. The numerical results and discussions are given 
in Sec. HI. The brief summary is presented finally in the 
last section. 



II. MODEL 

The interaction between the flux qubit and the LC 
resonator in the experiment 5!| is described by 
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where a\ a are the photon creation and annihilation op- 
erators in the basis of Fock states of the LC resonator, g 
is the flux qubit-cavity coupling constant. The RWA has 
not been employed here. The effective Hamiltonian for 
the flux qubit can be written as the standard one for a 
two-level system 
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where A and e are is the tunnehng coupling between the 
two persistent current states and the transition frequency 
of the flux qubit. e = /p($ — $o/2) with Ip the persis- 
tent current in the qubit loop, $ an externally applied 
magnetic flux, and $0 the flux quantum. In the above 
two equations, the Pauli matrix notations (7k{k ~ x,y, z) 
are used in the basis of the two persistent current states. 
Then the Hamiltonian for the whole system reads 

H = -{ea^ + Aa^) /2 + hiJr ia''a+ -\ 

+rig{a) + a)a, (3) 

where w^ is the cavity frequency. For convenience, we 
denote 



hujq ^ v/e^ + A2,tan6l = A/e 
Then the final Hamiltonian is ( ?i is set to unity) 

H = — -^ [cos{0)(Jz + sin(6')(Tj;] + w,. I a) a + - 

+g {a) + a) a^ (4) 

By introducing the new operators 

A — a + a,B — a — a,a = g/iOr 
we have 
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where e± = (w^ ± ojg cos6') /2. Note that the linear 
term for the original bosonic operator a^(a) is removed, 
and only the number operators A'^ A and B^B are left. 
Therefore the wavefunction can be expanded in terms of 
these new operators as 
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For A operator, we have 
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B operator has the same properties. Inserting Eqs. (6) 
and (7) into the Schro dinger equation, we have 
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where 
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In principle, all eigenvalues and eigenfunctions can be 
obtained in Eqs. (10) and (11). As before, to obtain the 
true exact results, the truncated number Ntr should be 
taken to infinity. Fortunately, it is not necessary. It is 
found that finite terms in state (7) are sufficient to give 
very accurate results with a relative errors less than 10~^ 
in the whole parameter space. We believe that we have 
exactly solved the JC model numerically. The numerical 
results are given in the next section. 

Solutions in the symmetry point- The spectrum in 
the symmetry point (e = ) is particularly interesting in 
experiments. In this case, Eq. (3) becomes 
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Associated with this Hamiltonian is a conserved parity 
H, such that [Hq, H] = 0, which is given by 



H — g*'ro-y/4gi7rJVg-iTr(T„/4 j^j _ ^f 



, TV = flTa + cr J2 + 1/2, (13) 



where N is the excitation number operator. H has two 
eigenvalues ±1, depending on whether the excitation 
number is even or odd. So the system has the corre- 
sponding even or odd parity. It is easily proven that the 
wavefunction (6) with even and odd parity is of the form 
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where ^+ (^-) is corresponding to wavefunction with 
even(odd) parity. Inserting Eq. (13) into Eq. (3) gives 
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The level transition is only allowed between the even and 
odd parity, i.e. E^ ^ i?. . The transition between 
the levels with the same parity is forbidden, E\ ' <f> 

-E- . The optical selection rules related to the parity 
have been discussed in the microwave- assisted transitions 
of superconducting quantum circuits [71 Il6|. 



III. RESULTS AND DISCUSSIONS 

Diaz et al diagonalize a restricted Hilbert space to a 
certain number of photon states (in the Fock basis) and 
obtained fitted parameters [5[ [l^l- The optimum fit of 
the experimental results within the present theoretical 
scheme gives Ip — 515nA,g/2Tr — 0.82GHz,uJr/2TT = 
8A3GHz,A/h = A.25GHz, very close to their values. 
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FIG. 1: (Color online) Theoretical spectrum of the flux qubit 
coupled to the LC resonator from numerical calculations. 
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FIG. 2: (Color online) (a) The energy levels EiT' and (b) 
the first 10 spectrums E„ — > ^'o ^^ ^^® symmetry point. 



The calculations in this paper are based on these param- 
eters, unless specified. 

We plot the numerical results for the spectrum for 
En — >■ Eo{ n = 1,2, and 3) in Fig. [TJ The experimental 
three spectral lines are just corresponding to the transi- 
tions between a few low energy levels, such as E3 — > Eo 
(upper), E2 — >■ Eq (middle), and Ei — >■ Eq (down). It is 
very interesting that our theoretical results for the spec- 
trum are in excellent agreement with the experimental 
ones in Fig. 3 of Ref. Q. Using these fitted parameters, 
the energy splitting on resonance (£'2 — Ei)/h obtained 
within the present approach is around O.OSTGiJz, just in 
the scope of the experimental observation. 

In Fig. 3 of Ref. [5|, a weakly visible spectrum line 
just below the middle spectrum line was attributed to 
the thermally excited qubit. We calculate the spectrum 
line for the transition E3 ^ Ei, as also list in Fig. [T]with 
a yellow line. Interestingly, it is just in the location ob- 
served experimentally shown in their Fig. 3. We believe 
that the state with Ei is just corresponding to the qubit 
excited thermally mentioned in Ref. [3] . 

We would like to mention here that the experimen- 
tally observed spectrum lines have been explicitly re- 
lated to the specified energy level transitions in the JC 
model without RWA. Then the comparison are easily per- 
formed. 

Next, we specially consider the case in the symmetry 
point . Fig. [2] (a) presents the energy levels from the nu- 
merically exact calculations. It was suggested in Ref. [5] 
that in the blue sideband spectral line [4] the minimum 
vanishes since the qubit is in the symmetry point where 
it produces no net flux and the transition is forbidden. 
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In the symmetry point, the transition from E. 
is forbidden due to the same parity, as shown in Fig. [5] 
This is the reason that the upper spectrum line around 
the symmetry point of Fig. 3 in Ref. [5| is almost invisi- 
ble. It is perhaps just the optical selection rules related 
to the parity makes the qubit to produce no net flux. 



As also indicated in Fig. [5J the other two transitions be- 
tween levels with the different parity are allowed, so the 
intensities in the middle and down spectrum lines in the 
symmetry point are nearly the same as in the whole spec- 
trum line. 



Fig. mb) shows the first 10 spectrums En 
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theoretically. In the experimental accessible detection, 
one can check the existence of these spectrums. 

We then turn to the Bloch-Siegert shift, which is just 
energy shift of the level transition with the consideration 
of the counter-rotating terms in the ultrastrong coupling 
regime. The Bloch-Siegert shift of the level transitions 
Ei — > Eq in the symmetry point are exhibited in Fig. 
[31 The smallest Bloch-Siegert shift is around 50MHz, 
nearly the same as that measured in the experiment i5i] . 
Note that some level transitions Ei -^ Eq are forbidden 
in the symmetry point due to the same parity, which 
are also presented here only for the estimation of the 
magnitude of the Bloch-Siegert shift in the corresponding 
spectrum. For the main spectrum Ei — > Eq, the Bloch- 
Siegert shift becomes larger as i increases, and its sign 
changes alternatively with either i, as shown in Fig. [3l^b). 

To show the effect of the qubit-cavity coupling strength 
on the Bloch-Siegert shift, we fix all parameters fitted 
from experiments except the coupling parameter g. The 
Bloch-Siegert shift of the level transitions Ei -> Eq in 
the symmetry point as a function of the effecting cou- 
pling constant a = g/i^r defined in Eq. (5) are plotted 
in Fig. m In the weak coupling regime, say a < 0.01, 
the Bloch-Siegert shift is so small ( less than 1 MHz) 
that it could not be distinguished from the spectrum 
line. When a > 0.1, the Bloch-Siegert shift increases 
considerably with the coupling constant, and can reach 
the regime of GHz. This observation is also of practical 
interest. Recently, the coupling could easily be further 
enhanced in the circuit QED |18l - l20l ] where g is compa- 
rable with ujr, i.e. a is in the order of magnitude of 
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symmetry point obtained within non-RWA (left panels) and 
RWA (right panels), (b) Bloch-Siegert shift for the different 
level transitions Ei — >■ £0 in the symmetry point. 
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observed to exceed 80MHz, the qubit line width at the 
symmetry point around AGHz, it is predicted that the 
Bloch-Siegert shift could be clearly resolved experimen- 
tally in this strong-coupling regime, like the Lamb shift 



IV. CONCLUSIONS 

In summary, by using extended bosonic coherent 
states, we solve the Jaynes-Cummings model without 
RWA exactly in the numerical sense. Within this tech- 
nique, we can reproduce excellently the spectrum mea- 
sured in a recent experiments on an LC resonator mag- 
netically coupled to a superconducting qubit [5J, which 
was demonstrated in the ultra-strong coupling regime. 
The Bloch-Siegert shift Ei — >■ Eq in the symmetry point 
is estimated to be 5QMHz, very close to the experimen- 
tal value. For the transition between the higher excited 
state i > 1 and the ground-state, the magnitude of the 
Bloch-Siegert shift monotonously increases, but the sign 
changes as (—1)'^*+-'^^. The considerable Bloch-Siegert 
shift in turn demonstrate that the counter-rotating terms 
should be considered. The effect of the qubit-cavity 
coupling strength on the Bloch-Siegert shift is also in- 
vestigated. It is predicted that the Bloch-Siegert shift 
can be distinguished experimentally for a > 0.15. The 
present technique are more suited for the stronger cou- 
pling regime, which experimental realizations may ap- 
pear in the near future [l8l-[2(j|. 
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